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I.  INTRODUCTION 


The  work  described  in  this  report  represents  one  phase  of  an 
experimental  study  of  the  electromagnetic  resonances  of  conducting 
bodies  with  attached  wires.  This  work  is  similar  to  a  previous  investi¬ 
gation  described  in  NASA  CR  169455  [Ref.  1];  the  major  difference  in  the 
present  case  is  the  use  of  smaller,  more  resistive  wires.  The  conducting 
bodies  included  two  cylinders  and  an  approximate  scale  model  of  an 
F-106B  aircraft.  The  results  from  the  cylinders  have  been  compared  with 
theoretical  calculations  to  check  the  accuracv  of  this  technique.  The 
results  from  the  aircraft  model  find  application  in  the  studv  ot 
lightning  strikes  to  airplanes.  The  wires  represent,  in  an  approximate 
sense,  the  lightning  channel.  Our  results  have  been  compared  with  those 
obtained  from  the  NASA  F-106B  during  direct  lightning  strikes.  The  need 
to  interpret  the  data  from  the  F-iOfeB  is  the  main  motivation  for  the 
work  reported  here. 

Section  II  of  the  report  describes  the  laboratory  technique 
emploved  to  investigate  the  resonances.  Short  pulses  of  current  were 
aoplied  to  the  bodv  under  test  through  one  of  the  attached  wires,  and 
f roe-field  electromagnetic  sensors  or  probes  were  used  to  measure  the 
B-dor  '  •)  B  M  r  )  and  D-dot  (  -ID/ It)  fields  as  a  function  of  time  near  tin- 
surface  of  the  bodv.  Two  wires  were  used,  one  for  current  entrv  and  the 
other  for  current  exit.  They  were  connected  axially  to  the  ends  of  tile 
cylinders  and  to  the  nose  and  tail  of  the  F-lOl.B  model,  with  the  curren; 
input  on  the  nose  wire.* 

A  curve-fitting  technique  known  as  Pronv  analysis  i  Re  t  s .  1  a 

w  a  -  used  to  stu  i  v  the  r-'snn  iirjs,  The  Pronv  code  was  run  on  the 
i  ,.ire  l  iut.i  *  ii-l  je'  s  of  'hi  j  e.-.  and  residues  were  ext  t  a<: "  ••• ; .  - 
•  ..  1  •  s  . i u  i  i  b~  ;  : i *  ■  >  r r  ■  r  - > i  •  i s  :•  he  natur  \  l  i  r  ■uu-'iit  z  * 

i.  -  in  i-tit”  .sv.i;  -ri.  Fourier  mai".,!..  -as  liso  useu  on  the  -lit.. 
i.  ■  .'iT.io;  lonroach  for  obtaining  information  >n  the  resonances . 


one  o 


The  Pronv  results  for  the  cylinders  are  given  in  Section  III.  Thev 
show  the  expected  weaker  damping  of  the  resonances  for  the  resistive- 
wire  case,  and  they  are  in  agreement  with  theoretical  calculations 
[Refs.  5,6]. 

The  Pronv  results  from  the  aircraft  model  are  given  in  Section  IV 
and  compared  with  results  obtained  on  the  NASA  F-106B  [Ref.  7].  The 
comparison  shows  that  the  use  of  resistive  wires  brings  the  resonances 
of  the  model  into  better  agreement  with  those  observed  in  flight. 

The  results  are  summarized  and  conclusions  drawn  in  Section  V'. 


II.  EXPERIMENTAL  SETUP  AND  MEASUREMENT  TECHNIQUE 


Experimental  Setup 

A  diagram  of  the  experimental  setup  is  shown  in  Figure  1 .  It 
consists  of  a  pulse  generator  (Tektronix,  Tvpe  109},  a  12-ft  bv  12-ft 
ground  plane,  a  sampling  oscilloscope  with  the  appropriate  plug-ins 
(Tektronix  Type  568  Oscilloscope,  Type  3S2  Sampling  Unit,  Tvpe  3T2 
Random  Sampling  Sweep  Plug-ins,  and  Type  S4  Sampling  Heads),  a  monitor 
oscilloscope  (Tektronix  Tvpe  7313),  some  in-house-built  buffer  amoli- 
fiers,  and  a  computer  with  floppv  disk  drive  for  the  digitizing  ant! 
recording  of  waveforms  (DEC  PDP  11/04,  Plessev  PM-XSil).  In  the  experi¬ 
ment,  the  object  undergoing  testing  is  either  a  cvlinder  or  in  F-106B 
model  located  10  ft*  above  the  ground  plane,  and  attached  to  the  rest  of 
the  exoeriment  with  wires  having  a  resistance  of  3.0  SI  to  the  foot  and  i 
diameter  of  0.01  in.** 

A  roughlv  rectangular  Dulse  with  a  1.2-ns-wide  base  and  a  rise  time 
and  a  fail  time  of  0.25  ns  each,  is  applied  at  the  ground  olane.  The 
oulse  propagates  up  the  lower  wire,  over  the  object  under  test,  on  up 
another  resistive  wire  attached  to  the  top  of  the  test  object,  and  from, 
this  wire  to  a  low  resistance  wire  attached  to  the  ground.  The  EM  fie  1.1 
near  the  test  object  is  measured  with  free- fie  id  sensors.  The 

required  for  a  portion  of  the  nulse  to  be  reflected  from  the  nose  .r:  ri  -■ 

test  obiect  down  to  the  ground  Diane  and  back  up  again  is  20  ns.  This 

gives  i  data  window  20  ns  wide  in  which  to  sample  the  waveform  be  fa:  »  : 

is  corrupted  bv  reflections. 

P a f a Acnuisit ion  System 

The  teem  i  si  t  ion  of  data  from  the  probes  is  done  iv  a  ,•  •.■-•s  m  r 
i  .  i  ai  Lv  modified  £or  this  r  tak  ■■  i .  it  t  -m.  •  ••  immar. :  ••  1- 

i  L  lose  one  and  an  analo"-to-di  ?it  ai  converter  wni<  h  ii  1 1  i  m  •  •  -  •_  v  .  "h 
converter  and  the  programmable  clock  are  born  standard  -r.;  t  tl  1  •• 


*  To  convert  f'“t  to  motors,  multinlv  bv  I.U43  COO  [>'):. 
'•r*  To  i-iiirvrt  in. dies  to  meters .  multiply  bv  2.  iiCd  CCO 


Figure  1.  Diagram  of  experimental  setup. 
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available  boards  which  plug  into  the  Q-bus  of  the  PDP  11  series 
puter.  The  A/D  converter  is  a  Data  Translation  DT  1712.  This  board  has 
a  single  12-bit  converter  with  8  differential  input  channels  multiplexed 
into  it.  The  differential  inputs  are  advantageous  because  the  probes 
have  two  outputs,  and  it  is  the  difference  between  the  outputs  that  is 
of  interest.  The  input  range  of  the  converter  is  from  -10  V  to  +10  V; 
this  makes  one  least  significant  bit  equal  to  4.88  mV.  The  maximum 
throughput  rate  is  35  kHz.  The  A/D  board  needs  3n  external  trigger  to 
mark  the  start  of  the  waveform  to  be  converted.  This  signal  is  sunolied 
bv  the  programmable  clock,  a  DEC  KW11-K. 

The  external  trigger  used  on  the  clock  board  is  the  "data  window" 
signal.  A  "data  window"  is  generated  bv  using  the  horizontal  sweep  of 
the  samoling  oscilloscope  to  saturate  a  simple  single  2112222  transistor 
amplifier.  The  output  of  the  saturated  amplifier  is  essential!"  an 
asvmmetric  squarevave  that  is  then  used  to  accomplish  three  different 
tasks.  The  first  is  triggering  the  programmable  clock  in  the  computer. 
The  second  task  is  "windowing"  the  data  that  the  computer  is  "se-ing." 
The  horizontal  sweep  for  the  monitor  scope  comes  from  the  sampling 
oscilloscope,  so  that  the  two  oscilloscopes  have  exac-lv  the  same  svoec 
r  f>>  and  the  traces  are  synchronized.  The  third  task  is  •»  a  t  i  n  g  a 
oounr-r  to  measu:  e  the  length  of  time-  it  takes  for  the  sampling  ns'iios- 
;  •*  i re-  •  T.-np  \  »w  —  -p  . 

The  display  on  ,  s  ire  :  in  '  osc  i  !  1  c  sour-  i the  r-sul:  of 
re  •  -c  I"  it-'  :>n  !  s-*  >.  Tli’is  *  u.i  1  ,.vv-r  rife  i.  nu<  h  si  -.v»»r  thin  r  ■ 


figure  (>.  Waveforms  recorded  for  the  small  cylinder. 
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Iv'jvcl  orm:;  recorded  lor  ihe  small  cylinder. 


III.  RESULTS  FROM  THE  CYLINDERS 


Measurements  of  the  natural  frequencies  of  the  cylinders  were  done 
first  for  comparison  to  previous  cylinder  work.  By  comparing  the 
present  results  with  those  of  Turner  [Ref.  1 j ,  the  degree  to  which  the 
resistive  wire  affected  the  experiment  was  found,  since  Turner  used  the 
same  cylinders  but  used  the  outer  shield  of  0.141-in  semirigid  coaxial 
cable  for  the  wires.  Comparisons  with  the  theoretical  calculations  bv 
Tesche  [Ref.  5]  and  Yang  [Ref.  b)  were  also  interesting.  The  work  bv 
Tesche  involved  the  natural  frequencies  of  isolated  cylinders,  while  the 
work  bv  Yang  dealt  with  the  effect  that  a  resistive  wire  attachment 
would  have  on  the  natural  frequencies  of  a  cylinder. 

Both  the  magnetic  field  (B-dot)  and  the  electric  field  (  D-dot  I  we.-e 
measured  at  the  center,  lengthwise,  of  the  cylinders.  As  shown  in 
Figure  5,  the  B-dot  probe  would  "see"  onlv  the  fundamental  frequent v  and 
its  odd  harmonics,  while  the  D-dot  probe  would  measure  onlv  the  o'.-n 
harmonics.  By  making  the  measurements  in  this  manner,  the  two  omf-s 
would  complement  each  other.  The  D-dot  probe  was  also  moved  to  a  second 
la  ation  on  both  of  the  cylinders  in  order  to  measure  the  odd  harmoni.s 
for  comparison  with  the  B-dot  results.  The  amount  of  agreement  of  the 
odd  harmonics  was  taken  as  a  measure  of  the  accuracy  of  the  tvchm  - . 
This  second  location  was  one-quarter  of  the  length  of  the  cylinder 
the  end. 

A  tvoicai  measured  B-dot  and  D-dot  response  for  the  small  c-linde.- 
is  shown  in  Figure  6.  uid  for  the  large  cylinder  in  Figure  7.  Tk  •  i 
pulse  that  was  used  is  ■>  iven  in  Figure  $.  Pronv  analvs  is  was  c  i-  : 
out  as  described  in  Reference  1.  The  onlv  special  processing  that  •; > 
waveforms  re.-.,  i  ved  before  bein'  ami":-.1,  bv  tin-  Pronv  ni  •  •  •.  • 
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Figure  4 


Photograph  of  the  model 


I-  i v> 1 1 r t-  i.  Sr ,1 1  ed  drawing  of  the  model  and  the  F-106B  aircraft 


by  Yang  [Ref.  6].  Some  direct  comparisons  with  their  calculated  poles 
are  in  Chapter  III. 

The  aircraft  model  is  an  approximate  model  of  an  F-106B  delta-wing 
aircraft.  The  model  was  constructed  in  the  following  manner.  The 
fuselage  was  made  of  an  aluminum  cylinder,  2  ft  long  with  a  4-in  dia¬ 
meter,  and  an  aluminum  cone,  1  ft  long  with  a  base  diameter  of  4  in, 
tapering  down  to  a  diameter  of  2  in.  The  tail  and  both  wings  were 
constructed  of  1  / 1 6-in-thick  brass  and  made  to  scale  with  the  rest  of 
the  model.  Thev  were  mechanicallv  attached  to  the  fuselage  with  screws, 
and  to  assure  a  good  electrical  connection,  copper  tape  was  also  used. 
The  overall  scale  of  the  model  was  1S.8:1.  A  comparison  between  rh^ 
model  and  the  actual  aircraft  is  shown  in  Figure  3,  and  Figure  4  shows 
the  model  in  the  experimental  setup. 


craft.  On  the  cvlinders,  various  positions  for  the  probes  are 
used . 

7.  Using  the  counter  (HP  5314a  Universal  Counter),  measure  the  time 
necessary  for  the  oscilloscope  to  complete  ten  sweeps  and  divide 
this  number  by  ten.  This  number  is  the  average  time  required  for 
the  oscilloscope  to  complete  a  sweep. 

After  these  preliminaries  are  completed,  type  in  the  command  R 
DATA.  The  machine  will  query  back  for  the  necessarv  information  before 
running.  The  name  of  the  output  file,  the  settings  on  the  sumo!  in'1 
o  s  c  i  1  I  o  s  c  ope  .  and  the  current  sweep  rate  of  the  s'si-m  will  be  inr  i.n-- 
tion  required  fop  the  program  to  proceed.  The  program  will  take  e. 
consecutive  sweeps,  deleting  the  first  sweep  and  keeping  the  last  ran. 
and  average  them  to  obtain  a  single  waveform .  During  the  comi.-ir 
sampling,  the  counter  should  be  left  on  in  order  to  measure  the  • 
necessary  for  the  eleven  sweeps  to  be  completed.  If.  due  to  a  malfunc¬ 
tion,  the  actual  rate  varies  from  the  rate  that  was  inserted  in  -h- 
program,  the  data  should  be  purged  from  the  records.  There  ia  si  wavs 
some  small  variation,  tvpicallv  0.10  to  0.25  percent,  which  is  accort- 
able . 

Sensors 

Two  different  sensors  were  used  to  make  all  the  measurement s  of  , he 
electromagnetic  field  on  both  the  model  and  the  two  cvl  ind-rs . 
were  a  D-dot  and  a  B-dot  probe.  The  B-dot  probe  is  a  model  MOL-cA .  . 

m  unt  far  t  tired  bv  CCiC.  having  a  bandwidth  of  at  l.-ast  1  . GHr  ; 

Tb.e  D-dot  probe  is  a  model  ACD-4At  R)  .  also  muf  tcfir<-*d  bv  2  Gad.  h 

i  bandwidth  of  i  r.  1  east  1.1  GHz  [  ">-f .  •»  i .  A  ••li-n  vr  aph  or  ‘  a  ■  ..ensor  • 

i  s  <j  i  ven  i  r.  F a  r  *  1 . 

first  evi  inuer  .  here  after  referred  to  as  sma  1  !  <:  v  !  •  n-n-r,  *..s 

long  and  bad  a  diameter  of  2  in.  The  second  •-  v  1  ind.'  r .  >•'  r:  a 

her  if  ter  as  r  he  lar’e  cvlinder.  was  also  3  ft  1  sn.> .  but  h  ad  a  ii  u-r.v 
or  t.  in.  The  ratio  of  d  i  ame  t  r- 1  o  -  1  en  g  t  h  of  the  small  c.  i  ind.-*;- 

.1 i  r  1  ”  tie*  as  rii.it  used  in  theoret  i  c  al  %ork  la  i  ■  :  ’>r.  *  ,r 


for  the  computer  to  digitize.  Taking  400  samples  in  this  period  of  time 
yields  an  actual  sampling  rate  of  206  samples / second . 

The  outputs  from  the  sampling  oscilloscope  have  an  impedance  of 
10  kft,  which  can  cause  a  problem  with  the  multiplexer.  If  the  cables  to 
the  multiplexer  have  too  much  capacitance,  there  is  an  undesirable 
"charge-up"  time.  There  are  two  wavs  to  correct  this  problem.  The 
first  is  to  use  short  cables,  but  there  is  a  limit  to  how  much  capaci¬ 
tance  can  be  removed  this  wav.  The  second  wav  is  to  lower  the  impedance 

feeding  into  the  multiplexer  bv  inserting  a  buffer  amplifier  in  the  iir.ir 
that  would  have  a  verv  high  input  impedance  and  a  verv  low  ofaur 
impedance.  The  high  input  impedance  of  the  amplifier  would  not  "lead 

down"  the  output  from  the  sampling  scope,  and  thus  eliminate  a  Dossibl*' 

source  of  distortion  of  the  waveform.  The  verv  low  outnut  imoedsnce  of 
the  amplifier  would  decrease  the  time  necess-iiv  to  charge  up  the 
capacitance  of  the  cables  and  the  assorted  strav  capacitances  in  th- 
eircuit.  The  second  wav  is  the  method  that  was  chosen. 

Experimental  Procedure 

Before  starting  the  data-taking  program,  the  following  procedure  is 

used . 

1.  Turn  on  all  the  equipment  (  except  the  pulser!  and  allow  it  at  r 

30  minutes  to  "warm  up,"  i.e.,  to  come  to  thermal  equilibrium. 

2.  Check  the  calibration  of  the  svarem  with  a  2  ns  standard  1:7  2. 

Time  Mark  Generator),  and  adjust  the  horizontal  sweeps  of  t 
oscil loscopes  if  necessarv. 

Turn  on  the  pulser  and  obtain  a  i r  of  signals  from  “he  or.j';,.. 

1.  Adjust  the  delav  in  the  B  channel  of  the  sampling  unit  plu-’-::;  L, 
that  the  two  signals  occur  s  i:i:ui  t  aneous  1  v . 
sir.  •  r.h“  DC  offset  ami  t  lie  !  1 . ;  •  ?•-  -  i  o:a  a  i  io.-tmeo'  s 

;  iaroiing  .uwoh  ,  adjust  -  a-  i  on  ,.f  “he  w  r’rfo  tin  it  the  ’  :v  • 

window  on  the  disoiav  of  t!a«:*  ::;or:  :  ■  .  r  usoj  ■  m.i  r*rt,ot  *.-  :  ,  ■ 

DC  level  . 

n.  Final  adjustments  of  the  position  of  the  probe  and  the  cables  fro:.: 

it  are  made  at  this  time.  The  probe  is  placed  at  a  position  r,e  i  > 

tr.e  model  that  correction  is  to  t  nos  it  ion  of  a  probe  on  r  he  a;r  - 


searching  for  a  minimum  in  the  frequency  spectrum  in  the  2-GHz  to 
2 . 5-GHz  band  and  making,  that  the  cutoff  frequency. 

The  Pronv  poles  have  been  normalized  in  the  following  wav.  The 
frequencies  (in  rad/s)  and  the  damping  rates  (in  Np/s)  were  multiplied 
bv  the  length  of  the  cvlinder,  L,  and  divided  bv  the  quantity  of  pi 
times  the  speed  of  light  in  vacuum.  See,  for  example,  the  labels  on  the 
axes  in  Figure  9.  The  normalized  frequency  of  the  first  resonance  thus 
has  a  value  near  1.0. 

The  results  of  the  Pron-'  analysis  on  both  of  the  fields  measured  at 
the  renter  of  the  small  cvlinder  are  given  in  Table  1 .  A  Pronv  order  f 
Id  with  a  sampling  rate  of  even'  sixth  point  was  used  on  the-  L-dot  dat.i . 
For  the  B-dot  data,  a  Pronv  order  of  13  was  also  used,  but  a  sampling 
rate  of  everv  eighth  point  was  used.  In  both  cases  the  program  was  set 
to  lake  ten  time  shifts.  Thus  the  real  poles  could  be  discerned  fro:,, 
the  pseudoooles,  created  by  the  Pronv  program,  on  the  basis  of  their 
stability.  Cnlv  the  poles  from  the  reconstructions  having  ax.  RMS  -rror 
Less  than  or  equal  to  6  percent  were  used  to  obtain  the  means  and  the 
standard  deviations  in  the  table. 

TABLE  1.  PRO. NY  RESLLTS  FOR  THE  SMALL  CYLINDER 
Pole  Number  Damning  i’reauenc v  ?:v:- 


0.231+0. DDL  0. "20+0.000  "-dot 

0.27 3+0 . 00 5  1 .d*  3+0.007  D- lor 

0. 304+0.0-,' 5  2.7  +  :*u.:oa  2-<i .  - 

0.12  >0.00  a  3.  u1,  7 +1-.00;  D-d,  t 


h-  re  woul-i  bv  some  overlap  of  fil->  'wi-j  tv  ,  ;.u  Oy  ;  .  s. •  r  '  . 

ills  of  this  ire  given  in  Table  2.  The  B  field  data  is  +  ne  s  vt- 
"  i  v“n  in  the  previous  table;  the  n«v  D-d. it  , :  at  a  h  is  i  i'r-tiv 


First 
Second 
Third 
Four  ti\ 


between  the  B-dot  and  D-dot  poles  is  generally  good  with  one  exception 
the  frequency  of  the  first  pole. 

TABLE  2.  A  COMPARISON  BETWEEN  THE  B-DOT  AND  THE  D-DOT  POLES 

OF  THE  SMALL  CYLINDER 


Pole  Number 

Damping 

Frequency 

Probe 

First 

-0.236+0.014 

0.851+0.004 

D-dot 

-0.231+0.004 

0.920+0.000 

B-dot 

Third 

-0.201+0.004 

2. 799+0.004 

D-iio  t 

-0. 30-+0.005 

2. 741+0.004 

3  - 1 1  rV 

The  different  v  dues  th.it  the  two  ne.isurejient  s  cave  for  r  i 
frequency  of  the  first  oo  le  was  disturbing.  Thera  were  two  po.«sihi 
sources  for  this  difference  in  the  two  waveforms.  The  first  sours.-  v>. 
that  the  probes  interacted  with  the  fields  of  the  cylinder  and  some iv 
either  raised  the  frequency  with  the  B-dot  probe .  or  lowered  th 
frequency  with  the  D-dot  probe.  The  second  possible  source  was  that  a 
accumulated  round-off  error  in  the  Pronv  prog  rim. 

To  examine  this  problem  further,  the  waveforms  we  re  "v  ace..  n.-d 
l.w-pass  filtering  program  that  w  as  1-si  -nv-J  :..j  ;u:t  ua^s  -  ,'i  ■ 

pole.  The  results  are  shown  in  Taoru  3.  The  j?  r-v.nent  ha  a  n .  . 
excellent.  From  this,  the  conclusion  is  drawn  that  diver  tei.ee  ot  : 

“  I  I*o  t  M  >  /  i  Va  r  j  VO  1  1  r*  r r  ?  r  . 

7API.Z  3 ,  A  DET  XT  LID  CCVVAIUZC':  OY  Tl!L  Fl?XT  POLE 


l  A  h* 


W .  i  £ 


the  cylinder;  then  the  D-dot  probe  was  moved  to  the  second  location, 
near  the  end  of  the  cvlinder,  and  the  D-field  was  recorded.  In  the 


analysis  on  the  first  set  of  data,  the  Prony  order  was  set  at  36  with  a 
sampling  rate  of  every  sixth  point  for  the  B-dot  waveform.  For  the 
D-dot  waveform,  the  Pronv  order  used  was  24  and  a  sampling  rate  of  everv 
sixth  point  was  used.  In  both  cases  the  acceptable  limit  on  the  recon¬ 
struction  error  was  set  at  6  percent.  The  results  are  given  in  Table  4. 

TABLE  4.  PRONY  RESULTS  FOR  THE  LARGE  CYLINDER 


Pole  N'umber 

Damping 

Freauencv 

Prnfc" 

F  ir  st 

-0.240+0.000 

0.327+0.005 

B-dot 

Second 

-0.200+0.000 

1 . 769+0.003 

D-dot 

Third 

-0. 330+0.024 

2. 722+0.06! 

B-dot 

Fourth 

-0. 347+0.005 

3 . 430+0.007 

D-dot 

For  the  analysis  of  the  D-field  waveform  measured  near  the  end  of 
the  cvlinder.  a  Pronv  order  of  30  and  a  sampling  rate  of  everv  sixth 
ooint  were  used.  The  results  are  <?  i  v«n  in  Table  3.  The  ig  re,-::  -nr 
between  the  odd  poles  obtained  from  the  two  probes  is  v«vv  good • 

TABLE  3.  A  '"CMPARL'CN  BETWEEN  THE  B-DOT  AND  THE  D-DOT  F'L.ES 

OF  THE  LARGE  CYLINDER 

Pole  Number  Damping  Freouen :v  P’"L 


-o.  2 wO +0.000  0.  >2  7+0.'!i 


comoarison  between  the  poles  generated  bv  the  two  >:v.  inde  i. 
o'-ovOj-d  Ln  tie*  7r  ,nh,  of  rL  ur.*  'L  In  this  u-inh.  "he  o  id  s  i  *s  ,,  ••  ;  ■■ 

:  r  i :•  ;.*■*  ■>  i  nd  D  f  i  elds. 


i 


O  • ) 


71  O 


difference  in  the  poles  is  the  lower  frequencies  of  the  larger  cvlin- 
der.  This  is  due  to  the  increased  capacitance  of  the  end  p'ates  of  the 
cylinder.  The  lowering  of  the  frequency  of  resonance  is  more  pronounced 
in  the  higher  modes.  A  second  effect  of  increasing  the  diameter  of  the 
cylinder  is  a  slightlv  stronger  damping  of  the  pole. 

A  comparison  with  the  results  of  Turner  I  Ref.  1|  for  both  cylinders 
is  given  in  Table  6  and  in  Figures  10  and  11.  In  Figure  10  the  small 
cylinder  results  are  compared,  and  in  Figure  11  the  large  cylinder 
results  are  compared.  The  wires  used  bv  Turner  veje  the  Conner  ■m'-.*; 
shield  of  0 . 1  i  1  -  i  n  semirigid  co.axia  1  cable;  r  h«- v  h  tve  j  ;nur  h  I  >  v.-  •- 
resistance  than  chose  u  >  e.j  ;  u  i  n-  present  •- !  u-  i  c:  w  >  •  i .  < .  j  !  ■ . 

diameter  - 

TAI'LF  0.  A  X.  Ml’A.cl  wlTd  P  kb .  !■'  '  ■'  V<  ' 

Turner  Virh  Resistive  Wir*  v?  *  »• 


Cmail 


0.  .  blvt-iO 


/♦■'ll  mot*4-*  so 

in  the 

L  h  • »  c 

1  i  f  f  e  r 

*  *  n  c  e  s 

.1  i  :  ih 

t  l  '  >  Tit 

f.he  *  * 

-.4  -.3  -.2  -.1 


<JL 

C  TT 


Figure  9.  A  comparison  between  results  from  the  two 
cylinders  with  resistive  wires. 


-.4  -.3  -.2  -.1 


a  L 

C  77 


igure  10.  A  comparison  between  the  small  cylinder 
results. 
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Figure  11.  A  comparison  between  the  large  cylinder 
results. 


wires  were  used  have  slightly,  though  consi stent lv ,  higher  damping 
rate.  The  effect  that  the  larger,  more  conductive  wires  have  on  the 
frequency  is  negligible. 

Theoretical  work  bv  Tesche  [Ref.  5]  covers  the  scattering  of  an 
electromagnetic  field  bv  an  isolated  cylinder  having  the  same  dimensions 
as  the  small  cvlinder  used  in  this  experiment.  Yang  (Ref.  6)  has 
performed  calculations  for  the  scattering  bv  a  cvlinder  which  has  a 
resistive  wire  attached.  The  ratio  of  the  dimensions  of  the  cvlinder 
i cvlinder  diame t e r / c vlinde r  length)  and  the  wire  attachment  'cylinder 
diameter  /  wire  diameter)  in  Van  * '  s  work  are  close  to  those  .if  :■ 
cvlinder  vised  in  this  exoeriment .  Van?  calculated  the  resonances  when 
the  resistivity  of  the  wire  was  2.51  Q/ ft  and  <  lien  it  was  251'J  L/fr. 
Yang's  model  could  calculate  oniv  the  odd  harmonics  of  the  cvl 
Table  7  gives  a  comparison  between  our  measured  results  for  the 
ind  the  results  the  two  computer  models  predicted. 

TABLE  7.  A  COMPARISON  WITH  SOME  THEORETICAL  WORE 
Pole  Number  Source  Damping  Normalised  Free. 


Vane  (2.51  11/ rt) 
Measured  (  •>  -.1;  ft  ! 


0.;52u 


i  \  •  .  t) 


;  a  u  :  *  ...  '  l  .  I  :  " 

Measured  '  : 

V  hi  ■  (  2  12  ■'  f- 


11“  ;  between  the  dampin'* 
T  V  •;  e  l.ne  re;;;;  l  f  3  C  1C. 


•’ll  ited  hv  Vine  r  f  no  two  d  i  r  f  •  r  v, " 
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This  section  on  the  resonances  of  cylinders  could  not  be  concluded 
without  a  comparison  of  the  results  of  the  Prony  analysis  with  the 
results  obtained  from  doing  a  fast  Fourier  transform  on  the  waveforms. 
Figures  13  and  L4  give  the  magnitudes  of  the  Fourier  spectra  of  the 
B-dot  and  the  D-dot  waveforms,  respectively.  The  resonances  are  revealed 
as  prominent  peaks  in  the  spectra.  The  locations  of  the  peaks  should, 
and  do,  agree  approximately  with  the  frequencies  of  the  poles  in  the 
tables.  For  example,  consider  the  poles,  -0.234  t  j0.880  and  -0.298 
i  j2.770,  listed  in  Table  7.  In  Figure  13(a),  the  peaks  corresponding 
to  these  poles  lie  at  about  0.158  GHz  and  0.452  HGz;  when  normalized, 
these  values  become  0.963  and  2.755,  demonstrating  the  approximate 
agreement.  Keep  in  mind  that  the  basic  frequency  resolution  of  the 
Fourier  transform  is  (20  ns)  or  0.05  GHz,  which  is  not  too  precise. 
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A  comparison  between  the  measured  results 
and  those  calculated  by  different  computer 
models . 
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I'uurict  s|ifi.ir.i  i>t  the  b-ilol  waveforms  Iron,  the  cylinders. 
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craft.  The  earlier  model  had  a  blunt  nose,  while  the  model  used  in  this 
experiment  utilized  a  tapered  one  to  achieve  a  more  exact  representation 
of  the  aircraft.  One  result  expected  from  tapering  the  nose  was  a 
slight  rise  in  the  damping  rates,  because  the  taper  of  the  nose  would 
act  as  a  transformer,  matching  the  impedance  of  the  rest  of  the  model  to 
the  wire.  Bv  matching  the  impedances,  the  reflection  coefficient  is 
lowered,  resulting  in  increased  damping,  of  the  waveform. 

As  can  be  seen  in  Table  12  and  in  Figure  22,  the  damping  rates  of 
all  but  the  first  pole  of  the  tapered  nose  model  were  subs  tan t ia  1  lv 
!  .'.we  red.  The  lack  of  significant  change  in  the  dampin’  rate  of  tV1 
first  pole  oav  have  been  rhe  result  of  the  offsetting  effects  of  rile 
resistive  wares  and  the  tapered  nose. 
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the  model  were  higher  than  those  from  the  airplane.  This  trend  was 
reversed  on  the  next  pole.  With  the  third  pole,  the  damping  rate  and 
the  frequency  of  the  pole  of  the  aircraft  was  higher  than  that  of  the 
model.  The  fourth  pole  of  the  aircraft  had  a  considerably  higher 
frequency  than  did  the  model,  and  a  moderatelv  higher  damping  rare.  The 
aircraft  did  not  have  a  pole  that  corresponded  to  the  fifth  pole 
generated  bv  the  model.  For  the  last  two  poles,  the  sixth  and  seventh, 
the  poles  from  the  model  had  a  slightly  higher  damping  rate  and  a  lower 
frequency  than  those  from  the  aircraft.  Overall,  the  approximate 
model  with  the  simple  wire  model  used  for  the  light  nine  ,rh  inne :  wo---. 
reasonably  well. 

TABLE  1  1  .  A  COMPAilTSO>;  BETWEEN'  THE  MODEL  AND  THE  AIRCRAFT  RE:' id. TO 
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Figure  20.  A  comparison  between  the  B-dot  and  the  D-dot 
results. 


calculating  the  damping  of  a  weak  pole  in  the  presence  of  a  strong,  pole. 
Onlv  the  frequencies  of  the  seventh  pole  will  be  compared,  and  the 
damping  rate  of  the  pole  from  the  D-dot  waveform  will  be  taken  as  the 
true  one.  The  fourth  pole  in  both  the  B-dot  and  the  D-dot  waveforms  is 
strong.  The  difference  in  the  poles  in  this  case  is  a  difference  of 
both  the  frequencies  and  damping  rate. 

The  two  sets  of  poles,  B-dot  and  D-dot,  are  displayed  in  the  graph 
in  Figure  20.  In  Table  10,  the  differences  in  the  two  sets  of  poles  are 
given.  In  this  table,  the  percent  difference  between  corresponding  pole 
pirts  is  calculated  as  the  difference  between  them  divided  bv  thei r 
average.  The  differences  are  seen  to  be  generally  quite  small . 

TABLE  10.  A  COMPARISON'  BETWEEN  THE  B-DOT  AND  THE  D-DOT  PRONY  RES!’ LIS 
Pole  Difference  in  the  Difference  in  the 


Numbe r 

Damping 

F requenc v 

First 

0.00 

0.400  MHz 

0 . 0 

5.326  R 

Second 

0.02 

0.000  MHz 

3.333 

0.603 

Thi  r»i 

o.o: 

0.270  MHz 

^ .  4  i )  :• 

i  .  -i  5  5 

Fa  u:r!i 

o.o  5 

2.  'CO  MHz 

21  .  277 

12.003  '■ 

2  i:;th 

0.00 

0.1  60  MHz 

0.0 

0 ,4^2 

f *v*nr  h 

Not  Compared 

! . 370  MHz 

: o  omr>  a  r  i  r.  ~  r  h-  Prone  r  •suit.,  t  rent  .  i  n  -  I  1  i  ’  h  t  v  .  v  ■'  ■  m 

:  Ref .  7  |  against  those  of  the  moiiel,  a  correlation  between  the  :ua  sets 
:  in  be  seen.  The  comparison  is  given  in  both  Table  11  and  F:  : u rj  2 1  . 
The  model  Poles  were  averaged  from  the  B-dot  and  D-dot  poles.  Lor  c  hi- 
first  two  po  t  h-*  d  inning  rate  and  the  r  reouonev  ef  rne  -o 
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poles.  The  program  examined  the  spectrum  of  the  unfiltered  waveform  and 
passed  only  the  frequencies  between  the  first  minimum  after  the  second 
pole  and  the  first  minimum  before  the  sixth  pole. 

For  the  waveform  that  contained  all  of  the  poles,  the  first, 
second,  sixth,  and  seventh  poles  were  obtained  along  with  a  combined 
third-fourth  pole.  The  combined  pole  was  not  used.  A  Pronv  order  of  24 
with  a  sampling,  rate  of  every  sixth  point  was  used  on  this  waveform.  The 
percent  of  error  (6  percent^  on  the  reconstruction  was  a  bit  higher  on 
the  Pronv  of  this  waveform,  but  the  poles  were  stable.  The  third  and 
fourth  poles  were  obtained  from  the  bandpassed  waveform.  The  Pr  re¬ 
order  used  with  this  waveform  was  20,  with  a  samoLing  r  ot  i-wv 
eighth  point.  The  upper  error  limit  for  the  reconstructions  used  ;o 
obtain  the  mean  and  standard  deviation  on  these  two  poles  was  1.3 
percent.  The  Pronv  results  are  given  in  Table  along  with  a  ••rarh  f 
the  fast  Fourier  program  in  Figure  19. 


TABLE  9.  PRONY  RESISTS  FOR  THE  MODEL  B-DOT  WAVEFORM 

Pole  Number  Damping  Scaled  Frequencv 

(  MHz  ) 
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figure  18.  Fourier  spectrum  of  the  U-dot  waveform  from  the  model. 


TABLE  8.  PROSY  RECl’LTS  TOR  THE  MODEL  D-DOT  WAVEFORM 


Pole  Number 

Damping 

Scaled  Frequency 
(  MHz  ) 

First 

-0.265+0.020 

7.308+0.082 

Second 

-0.240+0.005 

14.752+0.230 

Third 

-0.170+0.005 

18.420+0.045 

Fourth 

-0.206+0.005 

25.604+0.030 

Fifth 

-0. 340+0.007 

30.720+0.032 

S  ixt’n 

-0.106+0.012 

36.205+0. 106 

Seventh 

-0. i 56+0.023 

41).  70  !+0 . 7 

In  Figure  18,  the  graph  gives  the  frequence  spectrum  oht 
the  waveform  with  a  di’itai  fast  Fourier  routine, 
correlation  between  the  results  of  the  two  method.-.  «  ?ro;,v  ar. 

obtaining  the  frequencies  of  the  noles  in  the  D-dot  wave: nr 
Fourier  results,  five  of  the  poles  present  in  the  I’ronv  r 
distinct,  while  two  are  hidden. 
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cii.  i  i  c  a  tel  murcer.  When  toe  usual  filterin'  and  anal  vs  is  r 
i'.Vi  )■',  this  vi  v  ■’  f  ra  ,  the  first.  Sv  rend  .  >i’;r',.  >.id  :  u  t 

,;■!  .  •  1  '  v  7 h  t  tiiir  d  a  a  i  r  curt  a  r  -  -  -  v  -  -  v,  n  ' 

:  so  1  v  i  n  e  into  jot.ir  it-  poles.  Sc  art  in  at  i  Pr-nv  -.r  hr 

i  i*  i .  ;.i »■.  »  i  :r.;  •  t  i : i ■’  r  i  t  •  .> .  **.  ;  i»-*  L  i  r  “  .•  i'  .*> 1 « r>  n  2  •  10  1  ‘  .  •  •  •  ■  ' 

;-,r  n  ,  were  rec.ilvei  1 :  at-i’nle:  '.’it  :  !.-  t  h  i  r.j  tud  :  •>• 

iv,  • r  -  i  r  ■  P~  a  sinrle  ooj-.-.  fh<*  t  i-mu:  .  ■  :  +  h  .  •••:.!  :  -• 


i  .  .  .:>.t ..1 


-*  o  T  '.•*  *.<**n  1 . 1  1  -u>o[\  ^  i  :'t  t  twc  \ ;  ir  -  ••  •;  . 

Li.ni  -  .‘jn  C  h Pronv  •  ?.  r\.ier.  t  he  wav  r.h^n  j  tn  r-u; 

"wo  ool*?.1)  w.i:j  CO  sco  *r  i  *.  f  X  h  *  ::i  f  corn  l  h*  *  i**st  .  rlu*  v.iv  'fur 


.  i  h)  i*  o  *'  :*  i  nj 


i  nu  i  -i  r  »'i  ,3  b  irwii  .i.j.  j  l  i  ; 


;*  no.  r 


■  i  . 


: 

f 

1 

C-J 

r 

■M 

1 

' 

k 

6 

O 

Kn 

H”* 

. 

. 

I;.'.-" 

■  ■  ■ 

r.  h'“  i  ” 

1  ' 

,s 

l--  i 

'  '  " 

X-  V  ’  ♦  1 

•  -  i 

j 

■  i| 

*x  i.  1 

C  -  i  f 

CT  !  ‘  I 

•  ■: 

1 

t^TT 

r  -*  !  I 

-  - 

ll 

f  i"  «  1  f 

F  -  i.’  a 

|  1  ft 

i  - 

r  r 

€r~  !  1 

lO  5  1 

i  -  - 

-■ 

j 

*  ’  1 
..xJ.-  I 

|—  f 

A  ■ 

*•“  j 

i ' 

■ 

»  ■  :  •  • 

;  .  .  • 

r  •  *Jr  ■ 

*■ 

U  U— 

j 

i 

“■ 

r  ■■ 

K 

(- 

I- r 

! 

K  i" 

■ - 

- .  ■ 

U 

1  :  ■  ■ 

Figure  15.  TDK  of  the  model  in  Che  experiment. 


IV.  RESULTS  FROM  THE  F-106B  AIRCRAFT  MODEL 


Time  domain  ref lectometry  was  used  to  test  the  experimental  setup 
with  the  aircraft  model  in  place.  The  output  of  the  TDR  was  expected  to 
show  the  large-scale  structure  of  the  experiment:  the  junction  of  the 
50-52  cable  to  the  ground  plane/resistive  wire,  the  junction  of  the  model 
to  the  resistive  wire  at  both  ends  of  the  model,  and  the  junction  of  the 
resistive  wire  to  the  ordinary  wire  that  runs  back  to  the  ground. 
Because  of  the  dissipative  nature  of  the  resistive  wire,  the  fine  le*'.m 
of  the  model  was  expected  to  be  lost  in  TDR.  As  shown  in  Fi-m-e  !  : . 
these  expectations  proved  true. 

The  probes  were  positioned  near  the  model  so  that  t  iu-v  wou.a 
correspond  to  the  positions  of  the  equivalent  probe  on  the  aircraft.  Th  • 
D-dot.  probe  was  placed  at  the  underside  of  the  model  and  lust  dev-  - •>.  - 

tin  of  the  nose.  The  B-dot  probe  was  located  cn  the  tcpsid-  cf  re¬ 

model  just  above  the  seam  where  the  wing  loins  the  fusel  a ’e.  fd- 
t.j  '.act  ic  field  was  nonuniform  in  this  region,  and  the  dimer,  si ■  u\s  >•:'  •  ■ 

probe  were  of  the  same  order  as  the  gradient  of  the  field.  Becaus- 
this,  the  output  of  the  probe  corresponds  tc  the  average  field  :  n  a  i 
the  volume  of  the  probe.  Tvpical  waveforms  recorded  from  th  ■  B-cot  m  : 
h-  ior.  pruh.rs  are  riven  in  Figures  lu  and  17.  respective]-/. 

hi')  dam:  i  r.  ■  rate  of  each  ooi?  was  normalised  as  in  r  h..-  r  ••.  la 

•a  inter  .  out  th-  frequency  of  each  nole  was  scaled  downward  so  * 

direct  c  -.-i.i  t  r  i  son  could  ho  made  with  the  results  or  the  acru.i  1  »i:  -  :  "  . 
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or  ahtainini  the  mein  ind  st  md  ird  deviation,  on!  v  reconsi  ru.-r  ian  •• 
i .-rror  rata  less  than  a  percent  were  used.  The  results  ire  r v---n  ir. 


figure  14.  l  onrier  spectra  of  the  !)-dot  waveforms  from  the  cylinders. 


thev  would  have  a  rudder  and  elevons  like  the  real  aircraft. 

In  both  the  B-dot  and  the  D-dot  waveforms  there  was  a  zero- 
frequency  pole;  i.e.,  the  Pronv  program  extracted  a  pole  which  repre¬ 
sented  an  exponentially  decaying  term.  In  the  B-dot  waveform  the 
zero-frequency  pole  was  extremely  stable,  having  a  mean  value  of  0.05 
and  a  variance  of  zero.  In  the  D-dot  waveform,  however,  the  zero- 
frequency  pole  was  rather  unstable.  The  mean  of  its  damping  was  0.060, 
but  it  varied  from  0.03  to  0.09  and  had  a  standard  deviation  of  O.U22. 
This  z e r o-f requencv  pole  is  probablv  part  of  the  pulse  that  was  used  to 
excite  the  model.  On  both  of  the  cylinders,  the  Pronv  program  il..r 
extracted  zero-f requencv  poles  in  the  D-dot  and  B-dot  waveforms.  In  rh- 
cvlinder  results  there  were  two  zero-f requeue v  poles  rather  than  one . 


V.  CONCLUSIONS 


In  simplest  terms,  the  experiments  described  here  consisted  of 
producing  an  electrical  transient  disturbance  on  an  object  and  using 
frequency-spectrum  analysis  to  study  the  details  of  the  disturbance.  The 
object  was  an  airplane  model  with  attached  wires  in  the  laboratory, 
which  allowed  us  to  examine  some  aspects  of  the  real-life  problem  of  the 
F-106B  aircraft  in  a  lightning  strike.  Our  main  spectrum  analysis 
technique  (Prony  analysis)  gave  a  few  numbers  with  which  to  "character¬ 
ize"  the  object  under  test,  and  we  have  looked  in  particular  at  those 
numbers  which  tell  how  quickly  the  disturbance  must  damp  out.  The 
present  work  differs  from  that  done  earlier  [Ref.  I]  in  that  a  specific 
change  was  made  in  the  model — different  wires.  This  change  resulted  in 
improved  agreement  in  damping  between  the  model  and  a  particular  set  of 
lightning  data  for  the  real  F-106B;  so  that,  roughly  speaking,  we  may 
conclude  that  the  lightning  channel  was  more  like  the  wires  used  here 
than  like  the  previous  wires.  Our  basic  technique  of  excitation  of  an 
electrical  system  with  a  transient  input  and  the  characterization  of  its 
damping  properties  through  Prony  analysis  could,  of  course,  be  applied 
to  other,  nonelectrical  systems  as  well. 

More  specifically,  regarding  the  resonances  of  the  cylinders  we  can 
state  the  following  results  and  conclusions: 

1.  The  comparison  between  the  large-diameter  and  small-diameter 
cylinders  shows  slightly  stronger  damping  and  lower  frequencies  for  the 
large  one  (Fig.  9).  Isolated  cylinders  would  produce  the  same  result. 

2.  The  comparison  between  the  present  poles  and  those  of  Turner 
[Ref.  1]  shows  less  damping  in  the  present  case  (Fig.  10,11).  This  is 
expected  since  the  smaller,  more  resistive  wires  in  the  present  case 
result  in  less  current  conducted  away  from  the  cylinder. 

3.  The  comparison  with  the  poles  of  Yang  [Ref.  6]  shows  reasonable 
agreement,  and  the  comparison  with  Tesche  [Ref.  5]  shows  the  sort  of 
difference  expected:  less  damping  in  Tesche's  case  since  his  cylinders 
were  isolated  instead  of  wi re -connected  (Fig.  12).  These  comparisons 
give  confidence  in  the  basic  correctness  of  our  technique  for  determin- 
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ing  poles. 

From  the  F-106B  model  we  have  the  following  results  and  conclu¬ 
sions  : 

1.  The  comparison  between  the  poles  extracted  from  the  B-dot 
sensor  data  and  those  from  the  D-dot  sensor  shows  good  agreement  for 
most  of  the  poles.  More  specifically,  the  agreement  is  good,  within  10 
percent,  for  poles  1,  2,  3,  and  6;  it  is  fair,  21  percent,  for  pole  4; 
but  the  data  are  insufficient  for  a  good  comparison  on  5  and  7  (Fig.  20, 
Table  10).  Ideally,  the  poles  should  agree  exactly,  so  the  discrepancies 
that  are  observed  (which  are  usually  less  than  10%)  give  an  idea  of  the 
accuracy  of  the  poles  of  the  model. 

2.  The  comparison  between  the  poles  of  the  model  and  those  of  the 
actual  airplane  shows  rough  agreement,  the  damping  of  the  first  pole 
being  responsible  for  the  largest  discrepancy  (Fig.  21).  Poles  from 
only  one  lightning  event  on  the  airplane  were  used  for  this  comparison 
[Ref.  7].  Other  poles  have  been  obtained  from  airplane  data  and  may  be 
seen  plotted  in  Reference  1,  but  these  poles  are  less  reliable  because 
of  larger  quantization  errors  and  the  lack  of  simultaneous  B-dot  and 
D-dot  waveforms  for  corroboration  of  the  values. 

One  would  like  to  have  pole  sets  for  both  the  model  and  the 
in-flight  data  in  a  situation  where  the  attachment  points  were  known  to 
be  the  same.  Then,  differences  in  the  pole  sets  could  be  interpreted  as 
resulting  from  the  lightning  channel  having  an  impedance  either  higher 
or  lower  than  the  wires.  Thus,  something  would  be  learned  about  the 
channel  and  its  effect  on  the  resonances.  Because  the  attachment  points 
for  the  in-flight  lightning  event  used  here  are  not  known,  current 
conclusions  cannot  be  too  specific  regarding  the  channel.  However, 
rough  agreement  is  being  obtained  with  the  use  of  the  present  wires,  and 
some  of  the  existing  discrepancy  may  be  due  to  attachment  point  location 
variations  between  the  model  and  the  in-flight  situations. 

One  other  possible  source  of  discrepancies  between  the  poles  of  the 
model  and  the  airplane  is  the  shape  of  the  model:  it  is  not  an  exact 
scale  model  of  the  airplane.  Future  work  should  perhaps  be  aimed  at 
making  the  required  detail  improvements  to  the  shape. 
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3.  The  comparison  between  the  poles  of  the  present  model  and  those 
of  Turner  shows  the  effect  of  the  change  to  resistive  wires  and  a 
tapered  nose.  The  damping  of  all  the  poles  but  the  first  has  been 
reduced  (Fig.  22).  This  is  as  expected  and  is  the  same  effect  seen  in 
the  case  of  the  cylinders.  On  the  model,  this  reduced  damping  improves 
the  agreement  with  the  in-flight  results. 

4.  The  distribution  of  the  poles  in  the  complex  plane  is  different 
for  the  F-106B  model  than  for  the  cylinders.  Whereas  the  poles  of  each 
cylinder  lie  evenly  along  a  line  which  slopes  gently  to  the  left  (Fig. 
10,11),  the  poles  of  the  model  are  rather  scattered  and  show  a  tendency 
to  lie  farther  to  the  right  at  the  higher  frequencies  (Fig.  21).  This 
is  also  true  of  the  in-flight  poles  (Fig.  21). 

Some  comments  are  in  order  regarding  our  experiences  using  Prony 
analysis  on  laboratory  data,  in-flight  data,  and  computer  generated 
data.  For  computer  generated  waveforms  which  consist  of  several  damped 
sinusoids  without  noise  or  distortion,  the  Prony  code  works  very  well, 
extracting  the  correct  values  of  all  the  poles,  both  damping  and  fre¬ 
quency,  even  when  some  of  the  residues  are  very  weak  compared  to  others. 
In  some  cases,  the  frequencies  of  the  poles  can  also  be  picked  out  by 
inspection  of  the  Fourier  spectrum  of  the  waveform.  However,  in  many 
cases  the  spectrum  simply  does  not  reveal  the  weak  poles. 

In  the  Prony  analysis  of  waveforms  which  are  measured  rather  than 
computer  generated,  there  are  two  problems.  First,  the  Prony  code  often 
will  not  fit  the  waveform.  That  is,  the  RMS  error  between  the  actual 
waveform  and  the  one  generated  from  the  Prony  poles  anu  residues  is 
larger  than,  e.g.,  50  percent.  This  is  a  common  occurrence  in  the 
analysis  of  in-flight  waveforms.  When  it  happens,  the  poles  are  not 
used.  Second,  in  the  case  where  there  is  a  good  fit  (RMS  error  <  5/0,  a 
question  exists  as  to  whether  the  poles  are  really  the  true  natural 
frequencies  of  the  object  under  test,  or  whether  they  differ  from  these 
because  of  noise  or  distortion  in  the  measured  waveform.  One  example  of 
distortion  is  the  quantization  error  discussed  in  Reference  7,  which  was 
found  to  lead  to  incorrect  damping  rates  for  the  poles. 

To  gain  a  degree  of  confidence  in  the  natural  frequencies,  the 
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practice  has  been  to  analyze  simultaneous  B-dot  and  D-dot  waveforms  and 
make  a  comparison  of  the  resulting  poles.  If  they  agree  closely,  which 
often  happens  for  the  laboratory  data,  they  are  accepted  as  giving  the 
true  natural  frequencies  (including  those  of  the  input  waveform). 

One  method  tried  on  the  model  data  when  B-dot  and  D-dot  poles 
differ,  was  to  filter  out  some  of  the  poles  and  then  re-run  the  Prony 
code  on  the  filtered  waveform.  This  gives  the  code  a  simpler  waveform 
to  work  with  and,  as  described  in  Section  IV,  can  lead  to  better  agree¬ 
ment  between  B-dot  and  D-dot  poles. 

The  Prony  code  appears  better  suited  to  measured  waveforms,  which 
have  their  pole  frequencies  well  separated,  than  to  those  with  closely 
spaced  poles.  For  example,  the  pole  extraction  was  less  troublesome  for 
the  cylinder,  where  the  sensor  was  located  at  the  center  so  as  to  pick 
up  only  every  other  pole,  than  for  the  F-106B  model  with  its  many  poles. 

For  some  measurements,  the  correctness  of  the  natural  frequencies 
can  be  checked  in  another  way — by  comparison  with  theoretical  calcula¬ 
tions.  This  has  been  done  in  the  case  of  the  cylinders. 
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